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ABSTRACT 

We present near-infrared observations of the low-mass deeply-embedded Class 0/1 system 2MASS 
J17112318-2724315 (2M171123) in the B59 molecular cloud. Bright scattered light nebulosity is 
observed towards this source in the K 8 images, that seems to trace the edges of an outflow cavity. 
We report the detection of a low-luminosity protostar 2M17112255-27243448 (2M17112255) that lies 
-8" (-1000 AU) from 2M171123. This is a Class I system, as indicated by its 2-8 ztm slope and 
IRAC colors, with an estimated internal luminosity of ~0.3L Q . We estimate a mass of ~0. 12-0. 25 
Mq for this source, at an age of 0.1-lMyr. Also presented is detailed modeling of the 2M171123 
system. The best-fit parameters indicate a large envelope density of the order of — 10~ 13 gcm~ 3 , and 
an intermediate inclination between 53° and 59°. The observed if s -band variability for this system 
could be explained by slight variability in the mass infall rate between 2.5E-5 and 1.8E-5 M /yr. 
The protostar 2M171123 exhibits a rarely observed absorption feature near 11.3 /urn within its 10 
jum silicate band. We find a strong correlation between the strength in this 11.3 //m 'edge' and the 
iJ20-ice column density, indicating the origin of this feature in the thickness of the ice mantle over 
the silicate grains. 

Subject headings: circumstellar matter - stars: individual (2MASS J17112318-2724315) - - stars: 
individual (2M171 12255-27243448) - ISM: jets and outflows - ISM: abundances - 
stars: pre-main sequence 



1. INTRODUCTION 

The collapse of the cold molecular cloud from which 
stars are born requires the presence of rotationally sup- 
ported disks in order to conserve the angular momentum. 
On a timescale of 10 5 yrs, a central protostar is formed 
surrounded by an infalling envelope and an accreting disk 
(Hartmann 2000). Then begins the main accretion or T 
Tauri phase, during which the central protostar builds up 
its mass from the surrounding infalling envelope and ac- 
cretion disk. As gas accretes onto the star, an accretion- 
driven stellar wind develops that sweeps up infalling ma- 
terial from the surrounding cloud and drives a molecular 
outflow. The cavities carved out in the circumstellar en- 
velopes by such molecular outflows can be observed in 
scattered light nebulae. Several surveys based on near- 
infrared (NIR) and Spitzer/TRAC observations of outflow 
cavities in scattered light have enabled detailed study of 
the morphology of protostellar systems, such as the ex- 
tent of wind-created cavities, the large-scale geometry 
and inner structures of the infalling envelopes, as well as 
the size and structure of the disks (e.g., Padgett et al. 
1999; Burrows et al. 1996; Tobin et al. 2008). 

This paper discusses NIR observations and model- 
ing of the deeply embedded low-mass object 2MASS 
J17112317-2724315 (hereafter 2M171123) 'in the B59 
molecular cloud. B59 is an irregularly shaped nearby 
(d = 130^20 P c ! Lombardi et al. 2006) dark cloud and 
is the only known star-forming region in the Pipe Neb- 
ula. A recent Spitzer investigation of the dark cloud 
core under the "Cores to Disks Legacy Survey" (c2d) 
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resulted in doubling the previously known population, 
with 20 candidate low-mass young stars identified near 
the core, including two deeply embedded protostellar ob- 
jects (Brooke et al. 2007; hereafter B07). 2M171123 
was classified as a Class 0/1 system with a bolometric 
luminosity of 2.2±0.3L Q , based on its spectral energy 
distribution (SED) from NIR through millimeter wave- 
lengths. Extended structures possibly tracing the edges 
of an outflow cavity were observed towards this protostar 
in the IRAC 3.6 and 4.5 /zm images (B07). The apparent 
axis of the northeast extension seemed to be aligned with 
some of the CO outflow peaks earlier detected by Onishi 
et al. (1999; hereafter 099), suggesting 2M171123 to be 
the driving source of this outflow. Considering that no 
extended emission was observed in the 2MASS images, it 
was essential to obtain deeper NIR observations of this 
object to detect scattered light nebulae and conduct a 
detailed study of the outflow cavity. In $2] we present 
our NIR observations as well as some archival data uti- 
lized for this study, $3] discusses the morphology of this 
system as observed in the NIR images along with the 
detection of a new low-mass object ( 93.1]) . detailed mod- 
eling describing the structure of 2M171123 is presented 
in ij4j while $5] attempts to explain the 11.3 fim 'edge' 
observed in the 10 fim spectrum for 2M171123, thought 
to be indicative of the presence of crystalline silicates. 

2. OBSERVATIONS 

Observations in the J and K s bands of B59 cloud core 
were obtained at the CTIO Blanco 4 m telescope on 2008 
March 16, using the Infrared Side Port Imager (ISPI). 
This infrared camera has a field of view of 10.25' x 10.25', 
with a pixel scale of 0.3" pixel -1 . One pointing was de- 
fined at roughly a = 17^11™, <5 = -27°24', with ten 
dither positions obtained to minimize the effects of cos- 
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mic rays and bad CCD columns. The size of the dithers 
were between 20" and 30". The exposure times were 40s 
and 20s in the J and K s bands, respectively, with 3 num- 
ber of coadds obtained, thus resulting in a total time per 
image of 120s and 60s. Also presented here are H and K s 
band archival observations of Pipe nebula obtained at the 
VLT on 2002 July 28 and 30 [PID 69.C-0426(A)], using 
the Infrared Spectrometer and Array Camera (ISAAC). 
The SW imaging mode for ISAAC provides a field of view 
of 152" x 152", with a pixel scale of 0.148" pixel- 1 . The 
total integration times were 60s and 30s in the H and K s 
bands, respectively, with 10 dither positions obtained in 
each band. The dither sizes were between 2" and 10" 
in the if-band, and between 6" and 20" in the i^-band. 
The raw frames were bias-corrected and flat-fielded us- 
ing basic IR AF/ ccdred routines. For each object frame, 
a sky frame was constructed by median-combining all of 
the dithered observations excluding that object frame, 
without making any corrections for the offsets. The final 
sky-subtracted calibrated frames were first aligned to a 
common x and y position using the IRAF task imshift, 
and then averaged using the task imcombine. Aperture 
photometry was performed on the combined image for 
each band using the IRAF task apphot under the digiphot 
package. The data were calibrated using 2MASS mag- 
nitudes for bright sources in the same fields. The radial 
profiles for these bright sources were checked for any sat- 
uration effects. The magnitudes of the 2MASS sources 
that were used for calibration of the CTIO data ranged 
between J~11.5-11.9 mag and K s ~10.5-12.5 mag. For 
the VLT data, the magnitude range was between 13.6 
and 15.7 mag in the -ff-band, and between ~14 and 15 
mag in the -fT s -band. The calibration uncertainty is es- 
timated to be 0.08 mag and 0.04 mag in the ISAAG/H 
and K s bands, respectively, while it is 0.1 mag and 0.07 
mag in the CTIO/ J and K s bands. In Fig. [lj we show 
plots of our calibrated K s versus the 2MASS/-ftT s pho- 
tometry for the stars used for calibration. The dot- 
ted line in these plots indicates a linear least-square fit 
to the points, and we find a good correlation with the 
2MASS/ J FT S photometry for both the VLT and CTIO 
data. Stars that lie offset from the linear fit are repre- 
sented in Fig. [T] by open triangles, and have not been used 
as calibrators. These sources are 2MASS J17112508- 
2724425 (2MASS/iC s =11.57) in the CTIO data, and 
2MASS J17112111-2725004 (2MASS)/X S =14.99) and 
2MASS J17112749-2726029 (2MASS/# S =15.02) in the 
ISAAC observations. As mentioned, the radial profiles 
for these objects were checked for any saturation effects, 
and the linear fit shown in Fig. [T] indicates that satu- 
ration should not have affected the calibration for both 
the ISAAC and CTIO data. Table 1 lists the photom- 
etry thus obtained. 2M171123 could not be detected in 
the J-band, and an upper limit has been determined. 
We also list in Table 1 the 2MASS measurements for 
this target. In comparison with 2MASS/-?£T S , the CTIO 
measurement is ~1.2 magnitude brighter, while a less 
significant brightening of ~0.4 mag is observed in the 
ISAAC observations. [jH discusses the possible causes for 
the observed variability in the i^-band. fj5] compares 
the Spitzer/IRS spectra for 2M171123 with other pro- 
tostellar systems. The spectra were obtained from the 
Spitzer archives (PID 172, 179 and 20604). The spec- 
tra were extracted and calibrated using the Spitzer IRS 



Custom Extraction (SPICE) software provided by the 
Spitzer Science Center. For the purpose of modeling the 
2M171123 system, we have utilized Spitzer/IRAC and 
MIPS photometry from B07, as well as submillimeter 
and millimeter observations from Reipurth et al. (1996) 
and Wu et al. (2007). 

3. OBSERVED MORPHOLOGY 

099 detected 14 C 18 cores in the B59 region. One 
of the cores (named Core 1 in their paper) was found to 
have an emission intensity a factor of ^2 larger than all 
other cores. Two sources were found to be within 0.1 
pc of this C ls O core: 2M171123 and IRAS 17081-2721 
(2MASS J17111726-2725081; source 7 in B07). The CO 
outflow detected was also found to be closely associated 
with the C ls O peak (within 0.3 pc). The close prox- 
imity of 2M171123 and IRAS 17081-2721 thus suggested 
that one or both of them could be the driving sources 
of this outflow. 099 interpreted the outflow detected 
as a possible superposiion of two outflows; one outflow 
has a position angle of ~20° east of north, consisting of 
an elongated blue lobe ~0.7 pc and a compact red lobe 
within ~0.3 pc of the C ls O core. The other outflow has 
a position angle of ~ -40° and consists of a blue and 
a red lobe with spatial extents of 0.4-0.5 pc. Figure [2] 
shows the CTIO/AT s image, roughly centered around the 
2M171123 location. Also marked are the locations of the 
C 18 peak emission, the CO red and blue outflow peaks, 
as well as some of the sources from B07 discussed in this 
section. 

We observe a curved extended structure near 
2M171123 that extends to roughly 0.5 pc along the north- 
east direction. The extension is also clearly seen in the 
ISAAC/ K s observations (Fig. [3|). The apparent axis of 
this structure seems to be aligned with the first outflow 
discussed by 099, i.e. at a position angle of 20° east 
of north, and may trace the edges of this outflow cavity. 
This was first noted by B07 in their Spitzer/IRAC obser- 
vations, who gave a larger position angle of ~50° east of 
north. The boundaries of the eastern cavity as observed 
in the IRAC 3.6 and 4.5 /zm channels by B07 however 
seem more constraining for the position angle than the 
northeast nebula observed in our K s observations. 

Another extension is observed to the west of 2M171123 
in the CTIO/-ftf s image (Fig. [3^; marked by 'c'), with a 
spatial extent of about 0.007 pc. This faint emission 
could be due to the decreased sensitivity and worse see- 
ing (~0.8-l") in the CTIO data, as it is undetected in 
the more sensitive ISAAC/i^ image (Fig. [3b). It is also 
possible that the faint extension seen in the CT\0/K s is 
the northern boundary of the southwestern cavity, since 
this emission is at the expected location of that boundary 
(based on reflecting the southern boundary of the north- 
eastern cavity). The location of this structure is quite 
consistent with what is observed in the IRAC channels, 
although B07 have noted this extension to be confused 
by a superposed background source. This faint source 
(2M17112255-27243448; marked by 'b') is observed very 
clearly in the ISAAC /K s observations (Fig. [3Jd), and we 
discuss it in more detail in M3 - 1 1 

A more prominent curved structure is observed in 
the CTIO/# s -band towards the star 2MASS J17112153- 
2727417 (source 9 in B07) that seems to extend in a 
counter-clockwise direction (from bright to faint) first 



The Protostar 2MASS J17112318-2724315 



3 



Fig. 1. — 2MASS/-ftf a versus the calibrated _R" s -band photometry for the CTIO (left) and VLT (right) data. Dotted line shows a linear 
least-square fit to the data. Open triangles represent the stars not considered for calibration. 



TABLE 1 

Observations 



Star 


a (J2000) 


S (J2000) 




J 


H 




[3.6] 


[4.5] 


[5.8] 


[8] 


[24] 


2M171123 


17 11 23.18 


-27 24 31.5 


CTIO 


19.63 a 




13.87±0.11 


11.59±0.03 


9.30±0.03 


7.87±0.01 


6.50±0.01 


0.47±0.2 








ISAAC 




19.47±0.16 


15.46±0.07 


















2MASS 


18.80 a 


17.80 a 


15.08±0.14 












2M17112255 


17 11 22.55 


-27 24 34.38 


CTIO 
ISAAC 
2MASS 


20.62 a 


20.51±0.3 


15.98±0.14 
16.74±0.10 


14.2±0.10 


13.1±0.1 


12.7±0.2 


11.9±0.1 


3.66±0.2 



a Upper limit. 

towards the west, then northwest (Fig. [2]). This struc- 
ture may be associated with a third blue outflow peak 
(marked by an arrow) that lies off the image towards the 
southeast. However the apparent axis of this structure 
does not seem to align with this peak. The intensity of 
this blue lobe is weaker by a factor of ~2 compared to 
the red lobe, and lies outside the edge of the C ls O core 
(099). Source 9 lies —0.08 pc from this blue peak, but its 
closer proximity to the C ls O peak (~0.03 pc) suggests 
that this strong extension could be due to cloudshine 
emission. 

No clear extended emission is observed towards the 
red peak or the blue outflow peak west of 2M171123, 
suggesting that the second outflow at a position angle of 
-40° may be too weak in intensity to produce any strong 
scattered emission. Some strong nebulosity is observed 
towards sources 7, 13 and 14. A fainter structure of small 
spatial extent (—0.02 pc) was noted by B07 southeast of 
2M171123, that may be the southern or the southeast- 
ern boundary of the northeastern lobe. It seems that the 
southern boundary of the i4T s -band nebulosity is quite 
consistent with what has been observed by B07 in the 
IRAC bands. This structure however is not very promi- 
nent in the K s observations. 

3.1. Observed colors for 2M17112255 

In the CTIO and ISAAC images, we observe a faint 
source (a = 17 h ll m 22.55 s , 6 = -27°24'34.38"; here- 
after 2M17112255), -8" (-1000 AU) westward from 
2M171123. No 2MASS counterpart for this object could 
be found, and it was undetected in our J-band observa- 
tions. The NIR photometry along with a J-band upper 
limit is listed in Table 1. This is also a deeply-embedded 
object, with an A v ~ 19 - 30 estimated from the ex- 
tinction maps in B07. The c2d catalog lists a source, 



identified as SSTc2d J171122.5-272434, at the position of 
2M17112255. The 'object_type' for this source is given 
as 'YSOcred', i.e. a candidate YSO with a MIPS 24 
/im to IRAC flux ratio greater than 3. The c2d IRAC 
and MIPS flux densities are listed in Table 1. As dis- 
cussed in the c2d Legacy Project report 3 , to produce an 
IRAC-MIPS1 band-merged list, a 4.0" matching radius is 
utilized, while to merge this list with the MIPS2 (70 /im 
channel) source list, a larger radius of 8.0" is considered. 
The 24 /im photometry for 2M17112255 should thus have 
no contribution from 2M171123. There is no MIPS2 de- 
tection listed in the c2d catalog for 2M17112255. The 
MIPS2 detection quality for 2M171123 is listed as 'P\ in- 
dicating that there were multiple catalog detections near 
this source, which is expected given the 8" matching ra- 
dius. The 70 fim photometry is thus for the composite 
source. Fig 0k shows the SEDs for the two objects. The 
2 to 8 //m slope of the SED is measured to be 3.67 for 
2M171123, and 0.64 for 2M17112255. Following Wilking 
et al. (2001), we define the spectral index as a = d log 
(A-Fa)/ d log A, where a > 0.3 is a Class I source, 0.3 
to -0.3 is a 'Flat' source, -0.3 to -1.6 is a Class II source, 
and Class III sources have a < -1.6. Employing this 
criteria classifies both objects as Class I sources, as al- 
ready noted for 2M171123 by B07. Fig.gp compares the 
IRAC color-color diagram for the 20 low-mass candidate 
members in B59 with 2M17112255. The catalogued c2d 
photometry of 2M17112255 yields [5.8]-[8.0] = 0.82± 0.2, 
and [3.6]-[4.5]=1.06±0.14. In their study of IRAC colors 
of young stellar objects, Allen et al. (2004) have consid- 
ered the boundary between Class I and II objects to lie 
approximately at [3. 6]- [4.5] ~ 0.8 and [5.8]- [8] ~ 1.1. Ex- 
ceptions to this criteria are low-luminosity Class I objects 

3 available at |http: / /peggysue.as. utexas.edu/SIRTF/] 
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Fig. 2. — False color CTIO/i^s image, roughly centered around 2M171123 (marked by a rectangle). Also marked are the locations of 
C ls O peak (open square), red ouflow peak ('x'), and the blue peaks (plus signs). Also indicated are some of the sources from B07 discussed 
in [|3] The arrow points towards the third blue outflow peak that lies ^0.08 pc southeast from Source 9. The color intensity scale in units 
of MJy sr~ 1 is shown on the right. The image is 400" on a side, which corresponds to 0.25 pc at a distance of 130 pc. 



that would show bluer [5.8]-[8] colors at a given [3.6]-[4.5] 
color. Furthermore, Class I objects with low envelope 
densities of the order of ~ 10~ 14 (?cto~ 3 show bluer [3.6]- 
[4.5] colors at a given [5.8]-[8] color (Allen et al. 2004; 
Fig. 1). In comparison with the envelope models and 
the criteria for Class I objects discussed in Allen et al., 
the location of 2M17112255 in Fig.QJ) indicates that it is 
consistent with being a low- luminosity (~0.1Lq) Class I 
object with a low envelope density of ~ 3E — I4gcm~ 3 . 

For such embedded protostars, however, it is important 
to determine the internal luminosity, Lj n t, of the source, 
i.e., the total luminosity of the central protostar that 
excludes any emission from the external heating of the 
circumstellar envelope by the interstellar radiation field 
(e.g., Dunham et al. 2008). An estimate on Li nt could 
be obtained via 1-D or 2-D radiative transfer modeling of 
the observed SED. A detection at 70 (im though is crucial 
for such modeling, since the observed flux density at this 
wavelength is least affected by the source geometry and 
the presence or absence of a disk, and can better con- 
strain the resulting model-fits (Dunham et al. 2008). As 
mentioned above, the 70/im photometry listed in Table 
1 is for the composite source i.e. includes emission from 
both 2M171123 and 2M17112255. Our attempts at mod- 
eling the NIR to 24 ^m SED for 2M171 12255 resulted in 



largely degenerate model fits, mainly due to the sudden 
increase in the flux densities between 8 and 24 /xm, and 
the lack of data at longer wavelengths. We have instead 
considered the linear least-square relation between the 
observed flux at 24 /im and Li n t, obtained by Dunham 
et al. (2008) for a set of 11 embedded low- luminosity 
protostars. The relation is derived in the log-log space 
and is given as: 

log \F\ = 0.87±0.20 (log L int /L Q ) - 10.05 ±0.17, (1) 

where log XF\ is the flux at 24 /im in cgs units 
(ergs s~ 1 cm~ 2 ), normalized to 140 pc. Using this 
relation, L int for 2M171 12255 is estimated to be 
0.27tQ ^L Q . We note that while Dunham et al. have 
found a tight correlation between the 70/im flux and Li nt , 
the correlation at 24 /jm is weak (reduced x 2 value of 85, 
compared to 3 in the 70 fim case). Thus the estimate 
on Li n t for this source could best be considered as an 
upper limit. We consider the pre-main sequence evolu- 
tionary models by D'Antona & Mazzitelli (1994; 1998) 
to estimate the mass for 2M17112255 (Fig. gfc). The 
horizontal dotted line in this figure marks the estimated 
internal luminosity for this source. The age for B59 is 
not well-defined, and we consider the estimated range of 
0.5 - 1 Myr provided by B07. Within this range, the 



The Protostar 2MASS J17112318-2724315 



5 




Fig. 3.— False color images of 2M171123 in the CTIO/ifs-band (top, left), and the ISAAC/ifs-band (top, right). The labels a, b and 
c in the CTIO image indicate 2M171123, 2M17112255, and the faint bar-like extension, respectively. The color intensity scale in units of 
MJy sr~ 1 is shown on the right of each image. The images are 35" on a side, which corresponds to 0.02 pc at a distance of 130 pc. The 
bottom panel shows a A" s -band (2.16 fim) model image for 2M171123. The image is 25" on a side, corresponding to 0.016 pc at 130 pc. 
The overplotted contours correspond to intensity levels of 0.0006, 0.002, 0.005, 0.01, 0.015, 0.02 and 0.035 MJy sr' 1 . 



mass for 2M17112255 is estimated to be between ~0.18 
and 0.25 M©. If this object is as young as ~0.1 Myr, 
which is an age more consistent with it being a proto- 
stellar system, then its estimated mass could be as small 
as ^0.12 M Q . Considering the uncertainty on the inter- 
nal luminosity, logLi nt / Lq could be as small as ~ -0.8, 
placing 2M17112255 close to the sub-stellar boundary at 
an age of 0.1 Myr. This therefore could be the first can- 
didate brown dwarf in the B59 molecular cloud, though 
follow-up spectroscopy would be valuable in confirming 
its sub-stellar nature. 

We consider Poisson statistics to determine the prob- 
ability for 2M17112255 to be a background object. The 
surface density of young stars in the B59 core is found 
to be ~200 pc~ 2 (B07). This results in the expected 
number of stars within an 8" radius to be 0.0051, indi- 
cating a weak 0.5% probability of detecting a background 
source within 8" of a B59 member. Among the 20 can- 
didate low-mass stars in B59, there is then a 9 in 10 
chance of not finding a background star within this ra- 
dius. At a separation of ~8", 2M171 12255 thus has a 
weak 10% probability of being a background star or a 
chance alignment. Fig. [4jl compares the NIR colors for 
2M17112255 with 20 other candidate young stars in B59 
listed in B07. The observed intensity of 2M17112255 
in the low-extinction ISAAC/if s image may be consid- 
ered as an estimate of the intrinsic brightness of this 
object at this wavelength. In the H/ H-K color-color di- 
agram, 2M17112255 does not appear as a clear outskirt 
and follows the trend of redder H-K colors with fainter 
H. The zero-age main sequence (ZAMS) shown in the 
diagram was assembled from data presented in Bessell 



(1991) and Leggett et al. (1992). All of the 20 stars in- 
cluding 2M17112255 lie above the ZAMS, as expected for 
pre-main-sequence stars, with the more evolved Class II 
systems lying closer to the ZAMS than the less evolved 
ones categorized as 'Flat' and Class 0/1 or I sources by 
B07. The location of 2M17112255 thus suggests that it is 
unlikely to be a background field star. 2M17112255 and 
2M171123 could form a wide binary pair. It is interest- 
ing to note that the three known binary systems in B59, 
B59-1, B59-2 and LkHa 346, have projected separations 
>400 AU (Chelli et al. 1995; Reipurth & Zinnecker 1993; 
Koresko 2002), suggesting such wide systems to be more 
probable in this cloud core . However, a projected sepa- 
ration of ~1000 AU between 2M171123 and 2M17112255 
seems too large to form a stable gravitationally-bound 
system. The escape velocity for 2M17112255 from the 
primary is only ^0.6km/s, with a binding energy of 
~1E35 J (for a mass of 0.2M©). A weakly bound system 
such as this one may still be too young at an age of ~1 
Myr to be disrupted, even at such a large separation. 

4. MODELING 

We have used the two-dimensional radiative transfer 
code by Whitney et al. (2003) to model the 2M171123 
protostellar system. The main ingredients of the model 
are a rotationally flattened infalling envelope, bipolar 
cavities, and a flared accretion disk in hydrostatic equi- 
librium. For the circumstellar envelope, the angle- 
averaged density distribution varies roughly as p oc r -1 / 2 
for r << R c , and p oc r~ 3 / 2 for r >> R c . Here, R c is 

4 B59-1 is also known to have a tertiary companion at a projected 
separation of ~13 AU (Koresko 2002). 
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Fig. 4.— Top, left (a): SEDs for 2M171123 and 2M17112255. Top, right (b): IRAC color-color diagram for the B59 objects. Bottom, 
left (c): Evolutionary models from D'Antona & Mazzitelli (1998). The masses in Mq are marked. The isochrones (from top to bottom) are 
for 0.1 Myr, 0.5 Myr, 1 Myr, 10 Myr and 50 Myr. The horizontal dotted line indicates the internal luminosity for 2M17112255. Bottom, 
right (d): NIR color-magnitude diagram for the B59 objects. Bold line indicates the zero-age main sequence. 



the centrifugal radius and is set equal to the disk outer 
radius. It is the distance from the central protostar along 
the radial direction where material infalling from the en- 
velope lands onto the disk surface. A smaller value for 
R c would mean that high-density material would pile up 
closer to the protostar (e.g., Furlan et al. 2008). Thus 
increasing the value for this parameter would decrease 
the optical depth to the center of the envelope. 

The disk density is proportional to w~ a , where w is 
the radial coordinate in the disk midplane, and a is the 
radial density exponent. The disk scale height increases 
with radius, h = ha(za / 'i?*)^ , where ho is the scale height 
at i?„ and (3 is the flaring power. The disk density there- 
fore exhibits a Gaussian distribution about the disk mid- 
plane, i.e., falls off exponentially with scale height above 
and below this dense region, as well as in the radial 
direction towards larger radii. The disk extends from 
the dust destruction radius, R su b = R*(T su b/T*) 1 , to 
some outer disk radius, Rd,max- The dust sublimation 
temperature is adopted to be 1600 K. The disk is trun- 
cated sharply at its inner edge, and the region between 
the inner edge and the star is considered to be either de- 
void of any disk material or may contain some optically 
thin dust. 

Bipolar cavities are included in the models to account 
for the extended scattered emission observed at NIR 
wavelengths. The cavities extend from the center of the 
protostar to the outer radius of the envelope. We have 
adopted the curved cavity shape, the structure of which 
follows z = am 13 , where w = (x 2 + y 2 ) 1 ^ 2 . Here a is a 
constant determined by a relation between the envelope 



radius and the cavity opening angle, and b is the power 
of the polynomial defining the cavity shape. The shape 
parameter determines how quickly the cavity widens in 
the envelope. A small amount of dust is included in the 
cavity with constant density, nji 2 = 2 x 10 4 cm~ 3 (Whit- 
ney et al. 2003). 

Table 2 lists the best-fit model parameters for the 
2M171123 system. The stellar parameters used are il- 
lustrative: adequate for our purposes of modeling but 
not meant to be definitive. The values listed for these 
parameters are typical of T Tauri stars in Taurus, and 
were not varied during the modeling. We have used large 
grains in the dense disk midplane, with a size distribu- 
tion that decays exponentially for sizes larger than 50 /im 
up to 1 mm. The ISM-like grains with a max ~ 0.25 fim 
have been placed in the disk atmosphere and the outflow 
region. The grain model used in the envelope region is 
similar in size to the ISM-grains, except includes a layer 
of water ice on the grains that covers the outer 5% of the 
radius. 

Fig- El shows the best- fits obtained for our target. The 
main parameters that significantly effect the model SEDs 
are the inclination angle, the envelope density, and the 
centrifugal radius. As mentioned, decreasing the cen- 
trifugal radius implies that high-density material can ac- 
cumulate closer to the protostar. This results in more 
flux longward of ~20 /im and a deeper 10 /xm silicate 
absorption feature. The flux shortward of ~8 fim how- 
ever decreases since decreasing R c increases the optical 
depth to the center of the envelope, due to which lesser 
photons can escape through the cavity walls and the up- 
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TABLE 2 
Model Parameters 




Fig. 5. — Top, (a): The three different fits required to model the 
2M171123 system. Separate contributions from the disk, envelope 
and the star are shown in the middle panel for Fit#l, and in the 
bottom panel for Fit#2. 

per disk layers. Another important effect is the wave- 
length at which the peak in far-IR emission is obtained. 
Increasing R c shifts the region where infalling material 
piles up farther out to a larger radius, thus shifting the 
far-IR peak towards longer wavelengths. The MIPS/70 
pm observation is thus important in constraining the cen- 
trifugal and the outer disk radii, and a value of about 30 
AU for these parameters provides a good fit to the 70 
pm measurement for 2M171123. The effects of decreas- 
ing R c are similar to increasing the inclination angle, as 
both result in an increase in the amount of absorbing 
material in our line-of-sight. Therefore a smaller inclina- 
tion angle results in larger NIR fluxes, more emission at 
IRAC wavelengths, a shallower silicate absorption fea- 
ture, while a decrease in fluxes near ~100 /mi. Fig. [5] 
shows that while the 56° bin (Fit 1; covers inclinations 
between 53° and 59°) provides a good fit to the full depth 
of the silicate feature, the MIPS data points as well as 
the peak in the IRS spectrum, the fluxes remain small to 
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fit the CTlO/K s measurement. A better fit to this point 
(Fit 2) is obtained by slightly decreasing the envelope 
accretion rate from 2.5E-5 to 1.8E-5 M Q /yr, thus de- 
creasing the envelope density. Making the envelope less 
dense allows more photons to escape through the cavity 
regions, thus increasing the NIR flux. The silicate ab- 
sorption is also shallower for a less dense medium. The 
observed if s -band variability could thus be explained by 
a variable envelope accretion rate. This mass infall rate 
is related to the reference density pi, which is the den- 
sity the envelope would have at 1AU for the limit R c — 
(Kenyon et al. 1993): 



px = 5.3 x 10" 



M,, 



M, 



lO-SMQj/r- 1 / \IM Q 



-1/2 



(2) 



For the small range of infall rate for 2M171123, the ref- 
erence density is found to be between 2.4E-13 and 1.7E- 
13 gm/cm 3 . This is similar to the values obtained by 
Kenyon et al. (1993) for Taurus protostars, but is an 
order of magnitude larger than that reported by Furlan 
et al. (2008) in the same star-forming region. We note 
that 2M171123 is a deeply embedded source (B07) and 
lies close (^0.1 pc) to the C ls O peak that probes the 
densest, coldest envelope region, and so a larger enve- 
lope density is expected. 

Other parameters produce small variations to the 
model output. Among the disk parameters, more vis- 
cous energy is generated in the disk with an increased 
disk accretion rate, resulting in larger mid-IR fluxes. As 
the disk emission peaks near ^7 pm. for the viewing an- 
gle used, the parameter Mdisk is constrained by the 5-8 
/zm part of the spectrum. An Mdisk of 9.5xlO -7 M0/yr 
for 2M171123 is similar to the typical accretion rates ob- 
served among young protostars (~lx lO~ 8 M /yr; e.g., 
Whitney et al. 2003). The outer disk radius was set 
equal to the centrifugal radius. The inner disk radius is 
set by the star's magnetosphere. We set it equal to a few 
stellar radii, similar to the inner envelope radius. Since 
a cavity is included, the inner envelope radius will be 
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farther out from the central protostar. The outer enve- 
lope radius mainly effects the SED for A > 100 /im and 
can be constrained by the submillimcter and millimeter 
fluxes. However very large radii give too much optical 
depth to the center of the envelope, thus effecting the 
mid-IR fluxes. A value of —4000 AU provides a good fit. 

Fig. [3J; shows a 0.016x0.016 pc model image of 
2M171123, using the parameters from Fit 2. The model 
image shown has been obtained at an inclination of 60°, 
and was convolved with a Gaussian PSF of FWHM equal 
to that of 2M171123 in the CTIO/iT s observations. A 
foreground reddening of A v = 25.7 (B07) has been added, 
and the image has been rotated so as to match the ob- 
served northeast edge of the cavity. Making the system 
more inclined (i >80°) results in a bipolar cavity (e.g., 
Fig. 8 in Whitney et al. 2003), while decreasing it to 
~40° produces no outflow emission. Since reflected neb- 
ulosity is mainly observed in the northeast direction of 
2M171123 and the outflow is known not to be bipolar 
(099), a unipolar cavity produced at 50-60° inclinations 
is a better match to the observed morphology. 

The two main heating mechanisms in the protostellar 
system are accretion and protostellar irradiation. During 
the T Tauri accretion phase, the star's magnetosphere 
disrupts the accretion disk at its inner edge, causing ac- 
creting material to fall along the accretion columns at 
velocities of the order of —200 km/s onto the star (e.g., 
Hartmann 2000). The hot continuum emission seen in 
the optical and the ultraviolet wavelengths is produced 
by the accretion energy dissipated when the hot gas 
shocks at the stellar surface. The total accretion lumi- 
nosity, L acc , is therefore a sum of the viscous accretion 
energy generated by the disk, Ldisk, and the energy dis- 
sipated at the accretion shocks, L s h oc k- The total system 
luminosity is then 

Ltot — L s f ar + Z/s/jocfe + Ldisk- (3) 

Table 2 lists the luminosities from these three compo- 
nents. For a large inner disk radius as in the case of 
2M171123 {Rd.min ~ 36-R*), the energy dissipated at the 
hot spots will be the main heating source than viscous 
accretion within the disk. Almost 60% of the total en- 
ergy for 2M171123 is generated at these hot spots, that 
cover an area of 0.01 of the stellar surface, while only 
~1% is contributed by the disk. We note that a L to t of 
3.2L Q as obtained from the models is a factor of ~1.5 
larger than the bolometric luminosity of 2.2L Q reported 
by B07. Similar discrepancies have been noted by Furlan 
et al. (2008) in their modeling of Taurus protostars, and 
have been explained by the inclination angle of the sys- 
tem. Since the contribution from disk luminosity is min- 
imum for lowly inclined disks, model luminosity will be 
smaller than the bolometric luminosity. As the contri- 
bution from this component increases with higher incli- 
nations, models may overestimate the true bolometric 
luminosities (e.g., Whitney et al. 2003). 

None of the two models (Fits 1 and 2) at a small aper- 
ture of about 200 AU are able to reproduce the submil- 
limcter and millimeter measurements, and a much larger 
aperture of 5200 AU is required to fit these two points 
(Fit 3). The absence of any one good model that would 
fit these points along with the shorter-wavelength data 
indicates that the long-wavelength emission arises from 



a cold dust component that lies outside of the infall re- 
gion. The 350/im measurement corresponds to the flux 
integrated over a beam 40" in diameter (Wu et al. 2006), 
while the 1.3 mm observation is obtained from a 22" di- 
ameter aperture (Reipurth et al. 1996). We note that the 
millimeter observation plotted is actually of the source 
B59-MMSI that lies -15" from 2M171123, and therefore 
must include more emission from the surrounding dense 
molecular cloud than the protostellar system itself. The 
large aperture of 40" for the 350 /im measurement also 
indicates that cold dust from the surrounding dense cloud 
contributes a larger fraction to the observed emission at 
this wavelength. 

5. THICKNESS OF ICE MANTLE 

Recent surveys with Spitzer/IRS of a large number of 
protostellar systems have revealed a suite of absorption 
features, the most prominent being the ice features at 6.0 
and 6.8 firm, the 10 /im feature due to silicates, and the 
15 /im CO2 ice band (e.g., Boogert et al. 2008, Kessler- 
Silacci et al. 2005). At the low temperatures that exist in 
such dense molecular clouds and circumprotostellar en- 
velopes, atoms and molecules freeze out on dust grains, 
creating ice mantles that remain preserved in such cold 
environments. Thermal processes such as stellar irra- 
diation or shocks related to accretion activity or out- 
flows result in a thinning of these ice mantles, and alter 
the composition of the ice particles (e.g., Boogert et al. 
2008). 

Fig. [6]shows the optical depth, r, for the 10 /xm silicate 
and the 15 /jm CO2 absorption features for 2M171123. 
To measure the optical depth, the continuum in the 10 
/im silicate feature was approximated by a straight line 
that connects the 8 and 13 /im fluxes (e.g., Kessler-Silacci 
et al. 2005), while in the 15 /im feature it was approxi- 
mated with a straight line connecting the 14.65 and 16.3 
/tin fluxes (e.g., Quanz et al. 2007). The optical depth 
was then computed using 

F v /F c = e~ T , (4) 

where F v and F c are the observed and continuum fluxes, 
respectively. The silicate spectra for Class I objects are 
known to exhibit smooth, featureless absorption profiles 
with a strong, narrow peak at 9.6/im, indicative of amor- 
phous silicates (Kessler-Silacci et al. 2005). Stronger ab- 
sorption is observed in the more embedded protostars. 
A noticeable feature in the 2M171123 silicate spectrum 
is the significantly lower optical depth between —9.8 and 
11.3 /im. Such rare shallow absorption features have been 
reported previously for some low-mass embedded proto- 
stars, such as HH 46 IRS (IRAS 08242-5050), B5 IRS 
1 (IRAS 03445+3242) (Boogert et al. 2004) and L1489 
(IRAS 04016+2610) (Kessler-Silacci et al. 2005). This 
'shoulder' near 11.3 /im due to enhanced emission at this 
wavelength may be related to the presence of crystalline 
forsterites, commonly observed in older, more evolved 
Class II systems. The transition from amorphous to 
crystalline silicates results in a shift in the wavelength 
of maximum emission from 9.7 to 11.3 /im (e.g., Kessler- 
Silacci et al. 2005). Such transitions however occur on 
a timescale longer than the age of protostellar systems, 
which is why the presence of crystalline silicates in pro- 
tostars seems improbable. The main formation mech- 
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Fig. 6. — Optical depth in the 10 fim silicate feature (left) and 
the 15 /im C0 2 feature (right). From top to bottom: 2M171123, 
SST1711, HH 46, L1489 and B5 IRS 1. 

anism of crystalline grains in T Tauri disks is consid- 
ered to be thermal annealing in the warm inner disk re- 
gions, at temperatures >1000 K. Some low-temperature 
crystallization processes (e.g., Molster et al. 2001) may 
occur in these cold embedded sources, that may result 
in early crystallization of silicates present in the enve- 
lope. This would however suggest crystalline silicates to 
be ubiquitous among all protostars, and it is not clear 
why such processes may be active in just a few of these 
systems. The presence of a heating mechanism such as 
an outflow also could not have expedited silicate crystal- 
lization, since 2M171123, HH 46, B5 IRS 1 and L1489 
(IRAS 04016+2610) are all associated with an outflow 
source (Langer et al. 1996; Noriega-Crespo et al. 2004; 
Gomez et al. 1997), but display varying strengths in the 
11.3 /im feature (Fig. EJ). 

To understand this feature better, we have compared 
the 10 and 15 /im spectra for 2M171123 with a few other 
low- mass protostars, including B5 IRS 1, HH 46 and 
L1489. Among these, SSTc2d J171122.2-272602 (here- 
after SST1711) is another Class object identified in 
B59, and lies close to the CS molecular gas peak, as well 
as the extinction peak (B07). The 15 /im CO2 bending 
mode is an important tracer of ice composition, and the 
extent of ice processing that has occurred in the system 
due to various heating effects (e.g., Boogert et al. 2002). 



2M171123 and SST1711 show smooth, featureless C0 2 
profiles, similar to the spectrum of the background star 
CK 2 that probes the ice composition of the dark molec- 
ular cloud material, and is thus representative of unpro- 
cessed ice (Knez et al. 2005). HH 46 and L1489 show a 
double-peaked substructure at the base of the CO2 band 
that is caused by crystallization and effective segrega- 
tion of the H2O and CO2 ice species (e.g., Boogert et al. 
2004). The H 2 crystallization requires a temperature of 
at least 50K, while sublimation temperatures are higher 
(^90K; Boogert et al. 2002). The transition from amor- 
phous to crystalline ice thus occurs in this temperature 
range, suggesting that at least some fraction of the inner 
envelope for HH 46 and L1489 is warm enough (> 50K) 
to heat up the ice and cause the observed substructures 
(e.g., Boogert et al. 2004). The smooth C0 2 profiles 
for 2M171123 and SST1711 thus indicate temperatures 
below ~50K in their circumprotostellar envelopes and 
surrounding environment. 

An inverse trend is noticeable in Fig. [6] between the 
strength in the 11.3 /im shoulder and the optical depth 
in the CO2 band, with larger optical depths observed be- 
tween 9.8 and 11.3 /im for sources with deeper absorption 
at 15 /im. The origin of the 11.3 /im shoulder could be 
related to the thickness of the ice mantle coated on the 
silicate grains. As noted by Li et al. (2008), the inclu- 
sion of an ice mantle on silicate dust produces a shoulder 
at A > 11 /im, that may be interpreted as a crystalline 
silicate feature. A thicker ice mantle would then result 
in a more prominent shoulder at this wavelength. 

To confirm this hypothesis, we measured the strength 
in this feature, F pea k, and checked for any possible cor- 
relations with the water-ice column density. To mea- 
sure F pea k , we first approximated the continumm below 
the feature by a straight line that connects the observed 
fluxes at 9.8 and 11.3 /im, as shown in Fig. [7K (dot- 
ted line) for the case of 2M171123. These particular end 
points were selected since the shoulder is observed in this 
narrow wavelength range. The feature strength was then 
defined as: 

Fpeak = maximum [(F v - F c ,iq)/F Ci io + 1.0], (5) 

where f c ,io is the underlying continuum in the 9.8-11.3 
/im range, and (F v — i 7 'c,io)/^ r c,io are the normalized 
continuum-subtracted fluxes, plotted in Fig. [7b. Thus 
a larger value for F pea k would imply a greater contrast 
above the continuum, and therefore a more prominent 
shoulder would be observed. We repeated the same anal- 
ysis for the 15 /im CO2 absorption feature. In this case, 
the overlying continuum, F Cj i5, was approximated by 
connecting the observed fluxes at 14.65 and 16.3 /im, as 
shown in Fig. [7k. The normalized continuum-subtracted 
fluxes for the CO2 feature are plotted in Fig. [Tji. We de- 
fine a parameter, F15.2, as a measure of the contrast below 
the continuum -F C) i5, i.e. a smaller value for -F15.2 would 
imply deeper absorption in this ice feature. Using this 
method, we obtained F pea k and -F15.2 values for a set of 
low-mass protostellar systems, including 2M171123 and 
other objects shown in Fig. [S] These systems were se- 
lected based on the availability of the ice column density 
measurements, which have been obtained from Boogert 
et al. (2008). 

Fig. [8^ plots F pea k versus the water-ice column density, 
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Fig. 7. — Top, left (a): The observed 8-13 (im silicate spec- 
trum for 2M171123. Dotted line marks the continuum flux -F c .io 
approximated for the 9.8-11.3 (im shoulder. Top, right (b): The 
normalized continuum-subtracted fluxes in the 9.8-11.3 fim region. 
Bottom, left (c): The observed fluxes in the 15 fim CO2 ice-band 
for 2M171123. Dotted line marks the continuum flux F c 15 approx- 
imated in the 14.65-16.3 /im wavelength range. Bottom, right (d): 
The normalized continuum-subtracted fluxes in the 15 fim feature. 

N{H20). We note that both F pea k and F15.2 have val- 
ues in units of (F„ — F c .i ) / F c ^ + 1.0. An exponentially- 
increasing trend is observed between the two parameters, 
with F pea k increasing towards thicker ice mantles. F15.2 
shows an exponentially decreasing trend with N{H20) 
(Fig. [8p), as expected since a larger value for this pa- 
rameter indicates shallower absorption in the 15 fim ice 
band, and correlates with a smaller ice column density. 
In Fig. [5J;, F pea k is found to increase linearly with in- 
creasing depth in the ice feature, further confirming that 
objects with deeper ice absorption show greater strengths 
in the 11.3 /im shoulder. The lack of crystallization sig- 
natures in the CO2 absorption band for 2M171123 thus 
correlates well with the presence of a thick ice mantle and 
a more prominent 11.3 /im edge. Thus protostars that 
are not deeply embedded or have heated their envelopes 
sufficiently to melt away a large fraction of the ice coat- 
ings exhibit a weaker strength in the 11.3 /im shoulder, 
confirming the origin of this feature in the thickness of 
ice coating over the grains. We note that some crys- 
talline silicates may be present in the circumprotostellar 
envelopes, as ISM-like dust is known to show a 3-5% 
crystalline mass fraction (e.g., Li et al. 2008). Neverthe- 
less, it is the thickness of the ice coating over the silicates 
grains that seems to effect the observed 11.3 fim shoulder 
strongly. 

6. CONCLUSIONS 

We present NIR observations of the low-mass proto- 
star 2M171123 in the B59 molecular cloud. Scattered 
light nebulosity possibly tracing an outflow cavity are 
observed towards this source. Modeling this system in- 
dicates slight variability in the mass infall rate between 
2.5E-5 and 1.8E-5 M Q /yr, that could explain the ob- 
served i^-band variability. The system is viewed at an 
intermediate inclination of 56° to the line-of-sight. We 
report the detection of a low-luminosity Class I object 
2M17112255, that lies -8" from 2M171123. We estimate 
a mass of ~0. 12-0. 25 Mq for this source, at an age of 0.1- 
lMyr. 2M171123 shows a smooth, featureless ice profile 
with no signs of ice processing, indicating a circumproto- 
stellar envelope at a temperature below ~50 K. We find 
a strong correlation between the thickness of ice coating 




F peak 

Fig. 8. — Top, (a): F pea k versus the /^O-ice column density. 
Middle, (b): -F15.2 versus the ice column density. Bottom, (c): 
Fpeak versus F15.2. The thicker the ice mantle over silicate grains, 
the more prominent is the 11.3 /im feature. 

over the silicate grains, and the strength in the rarely ob- 
served 11.3 fim absorption feature. The enhanced 11.3 
/im emission observed towards 2M171123 must be due to 
a thicker ice mantle over the silicate grains. 
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